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The  complex  nature  of  lithium  ion  batteries  makes  the  tracking  of  individual  physical  processes  difficult. 
However,  having  a  good  knowledge  thereof  is  indispensable  for  a  continuous  improvement  of  cell 
performance  and  lifetime. 

In  this  study,  the  impedance  response  of  a  commercial  18650  cell  was  investigated  within  a  wide 
frequency  range  (100  kHz-2  pHz)  by  combining  electrochemical  impedance  spectroscopy  and  time 
domain  measurements.  In  parallel,  the  original  18650  cell  was  opened.  Thus,  the  electrode  materials 
could  be  tested  in  an  experimental  test  set-up  using  a  reference  electrode,  which  allowed  the  separa¬ 
tion  of  anode  and  cathode  processes.  The  different  kinds  of  impedance  data  sets  were  de-convoluted  in 
the  space  of  relaxation  times,  enabling  a  refined  separation  of  physical  processes  like  charge  transfer  or 
solid  state  diffusion.  This  multistep  approach,  which  is  not  tied  to  a  particular  electrochemistry,  allowed 
the  identification  of  the  dominating  physical  processes  being  hidden  in  the  impedance  spectrum  of  the 
original  18650  cell. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Lithium  ion  batteries  receive  growing  attention  because  of  their 
potentially  high  energy  and  power  density,  which  are  required  for 
“Electric  Vehicles”  (EV)  and  “Hybrid  Electric  Vehicles"  (HEV). 

However,  material  chemistry  and  fabrication  technology  of 
previous  cell  generations  must  be  improved  considerably  for 
application  in  electric  mobility.  This  task  requires  characterization 
methods  which  are  capable  to  assign  performance  characteristics  to 
distinctive  cell  components. 


Abbreviations:  DRT,  distribution  of  relaxation  times;  EC,  ethylene  carbonate; 
ECM,  equivalent  circuit  model;  EIS,  electrochemical  impedance  spectroscopy;  EMC, 
ethyl-methyl-carbonate;  EV,  electric  vehicle;  FIB,  focused  ion  beam;  HEV,  hybrid 
electric  vehicle;  KK,  Kramers  Kronig;  SEI,  solid  electrolyte  interphase;  SEM,  scan¬ 
ning  electron  microscope;  SOC,  state  of  charge;  TDM,  time  domain  measurement; 
VC,  vinylene  carbonate. 
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The  relevant  loss  processes  in  lithium  ion  cells  as  well  as  in 
alternative  cell  chemistries  occur  in  a  wide  frequency  range  from 
pHz  to  MHz  frequencies  [1],  For  instance,  the  solid  state  diffusion 
in  the  electrodes’  active  material  is  a  very  slow  loss  process 
occurring  at  frequencies  between  pHz  and  Hz  frequencies.  These 
frequencies  can  be  investigated  adequately  by  pulse  measure¬ 
ments  which  belong  to  the  group  of  time  domain  measurements 
(TDM).  Hereby,  the  system  response  of  lithium  ion  cells  is  evalu¬ 
ated  by  Fourier  transform  in  order  to  derive  the  low  frequency 
impedance  spectrum  from  current  pulses  [2],  For  very  high 
frequencies  above  5  kHz,  contact  resistances  between  electrode 
and  current  collector  or  between  the  active  material  particles 
contribute  to  the  impedance  response.  These  frequencies  are 
primarily  analyzed  by  electrochemical  impedance  spectroscopy 
(EIS)  [3],  which  is  also  used  to  analyze  the  medium  frequencies.  For 
these  frequencies  between  1  Hz  and  5  kHz,  losses  as  charge 
transfer  or  solid  electrolyte  interphase  (SEI)  are  dominating  the 
impedance  contributions. 
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A  combination  of  both  methods,  EIS  and  TDM,  was  recently 
introduced  for  the  characterization  of  lithium  ion  cells  by  our  group 
[2]  and  is  applied  in  this  contribution. 

However,  the  individual  time  constants  of  different  loss 
processes  are  often  close  to  each  other  and  therefore  prevent  the 
separation  of  more  than  two  polarization  processes  by  analyzing 
the  Nyquist  plot. 

It  was  already  shown  in  previous  papers  for  solid  oxide  fuel  cells 
[4—6],  that  it  is  advantageous  to  de-convolute  the  impedance  data 
in  the  space  of  relaxation  times.  This  approach,  called  DRT  (distri¬ 
bution  of  relaxation  times)  method,  enables  a  refined  separation  of 
physical  processes  like  charge  transfer  or  diffusion  losses.  Thereby, 
every  physical  process  is  represented  as  local  maximum  in 
a  continuous  distribution  function.  The  advantages  of  the  DRT 
method  are  not  tied  to  a  particular  electrochemistry. 

Recently,  it  was  applied  for  the  first  time  to  lithium  ion  exper¬ 
imental  cells  [7-9],  Applying  this  procedure,  cathode  and  anode 
polarization  processes  with  close  time  constants  were  separated 
individually  and  a  physically  motivated  equivalent  circuit  model 
(ECM)  was  established  without  any  a  priori  settings  for  the  elec¬ 
trochemical  system.  The  loss  processes  for  the  investigated  lab 
scale  LiFePCU  cathodes  have  been  pegged  as:  (i)  solid  state  diffusion 
in  the  cathode,  (ii)  charge  transfer  resistance  between  cathode  and 
electrolyte,  and  (iii)  contact  resistance  between  cathode  and 
current  collector.  This  allowed  for  the  determination  of  tempera¬ 
ture  and  SOC-dependency  for  each  individual  loss  processes.  The 
contact  resistance  shows  no  distinctive  parameter  dependency 
whereas  the  other  losses  show  pronounced  parameter  dependency 
on  both  parameters. 

An  application  of  this  method  on  commercial  cells  and  its 
electrodes  is  desirable  in  order  to  understand  and  afterward 
improve  the  current  cell  generation. 

Today,  the  analysis  of  commercial  cells  is  usually  based  on 
electrochemical  impedance  spectroscopy  at  varying  operation 
parameters  as  state  of  charge  (SOC)  or  temperature.  For  the  analysis 
of  parameter  dependency  of  loss  processes,  the  Nyquist  plots  of 
impedance  curves  are  complemented  by  the  analysis  of  Bode  plots, 
as  the  frequency  information  is  not  provided  in  the  Nyquist  plot 
visualization  [10,11],  The  loss  processes  are  usually  interpreted  by 
transferring  the  experience  of  former  investigations  or  publications 
to  new  results  [1,10,12],  According  to  these  identified  loss  processes, 
which  are  usually  not  more  than  two  polarization  losses,  equivalent 
circuits  are  proposed  in  order  to  evaluate  their  parameter  depen¬ 
dencies  [10,12],  The  standard  frequency  range  of  these  impedance 
evaluations  is  limited  to  10  kHz  for  high  and  10  mHz  for  low 
frequencies. 

A  more  detailed  analysis  of  number  and  physical  origin  of  loss 
processes  would  be  profitable  in  order  to  understand  the  reasons 
for  cell  performance  and  aging  mechanisms  and  to  avoid  the 
cumulated  analysis  of  several  overlapping  loss  processes. 

There  are  a  few  publications,  using  a  reference  electrode  which 
can  be  inserted  into  18650  cells  in  order  to  separate  the  major 
contributions  of  anode  and  cathode  [13,14],  Another  way  of  analysis 
technique  which  is  introduced  in  literature  is  to  open  the 
commercial  cell  and  to  analyze  the  single  electrodes  in  experi¬ 
mental  cells  by  reference  electrode  [15—17]  or  symmetrical  exper¬ 
imental  cells  [18].  Hereby,  an  analysis  of  single  electrode 
contributions  to  the  increase  of  impedance  polarization  or  to  the 
over  voltage  is  enabled  [17,18],  This  analysis  is  usually  performed 
for  frequencies  between  10  kHz  and  10  mHz.  The  evaluation  of 
smaller  frequencies  is  overlapped  by  the  intercalation  capacity 
which  allows  only  a  qualitative  evaluation  of  the  real  part  values. 
Furthermore,  the  frequency  range  is  limited  as  the  measurement 
time  increases  disproportionately  for  decreasing  frequencies.  The 
analysis  of  Nyquist  plots  allows  thereby  the  comparison  of  the 


entire  polarization  of  anode  and  cathode,  but  not  the  separate 
comparison  of  electrodes  losses  for  individual  frequency  ranges. 
The  information  about  time  constants  of  individual  loss  processes 
cannot  be  visualized. 

In  this  study,  the  impedance  response  of  a  commercial  18650 
cell  is  characterized  by  combining  EIS  and  TDM  as  shown  sche¬ 
matically  in  Fig.  1.  This  delivers  information  over  a  wide  frequency 
range  of  11  decades  (2  pHz-100  kHz).  Subsequently,  the  DRT 
method  de-convolutes  the  impedance  spectra  in  the  space  of 
relaxation  times  in  order  to  analyze  the  relevant  frequency  ranges 
separately.  Our  aim  is  to  identify  the  most  prominent  loss  processes 
which  limit  performance  or  lifetime  of  a  commercial  18650  cell. 
Furthermore,  the  origin  of  these  processes  shall  be  assigned  to 
anode,  cathode,  or  electrolyte,  or  interfaces  thereof. 

For  this  purpose,  a  widely  used  type  of  18650  cell  (LiFePO/t/ 
graphite)  is  first  measured,  then  opened  in  a  glovebox  and  the 
extracted  electrodes  are  re-analyzed  in  experimental  cell  set-ups. 
The  contributions  of  cathode  and  anode  are  separated  using 
a  reference  electrode.  Furthermore,  the  evaluable  frequency  range 
is  expanded,  as  impedance  spectra  of  experimental  cells  are  less 
spoiled  by  inductive  contributions  compared  to  18650  cells.  This 
allows  for  a  comparative  study  of  measured  impedance  spectra  and 
the  calculated  continuous  distribution  of  relaxation  times  (DRT)  of 
experimental  cells  with  the  original  18650  cell. 

New  insights  into  the  internal  loss  processes  can  be  given  as  DRT 
visualization  directly  provides  information  about  the  characteristic 
frequencies  of  loss  processes.  This  allows  for  an  independent 
evaluation  of  various  frequency  ranges  and  a  detailed  comparison 
of  experimental  and  commercial  cell  measurements  via  polariza¬ 
tion  and  characteristic  frequency  of  the  occurring  loss  processes. 

2.  Experimental 

2.1.  18650  Cell  and  its  components 

In  this  study  we  investigate  a  widely  used  type  of  18650  cell 
with  a  nominal  capacity  of  1.1  Ah.  The  cell  is  composed  of  LiFeP04 
and  carbon  black  as  cathode  and  graphite  as  anode  material  (Fig.  2). 
The  cathode  microstructure  shows  interesting  features  in  terms  of 
agglomerate  size  as  well  as  pore  size  distribution,  which  was 
extensively  studied  by  focused  ion  beam  (FIB)  tomography  in 
reference  [19],  The  subsequent  characterization  in  experimental 
cells  necessitates  opening  of  the  fully  discharged  18650  cell  inside 
a  glovebox  and  extracting  anode  and  cathode  sheets  thereof.  The 
following  procedure  is  different  for  anodes  and  cathodes. 

A  straightforward  extraction  is  possible  for  the  anode,  as  some 
parts  of  the  graphite/copper  sheet  are  coated  only  on  one  side  with 
anode  material:  1.  Punch  this  part  of  the  anode  sheet  with  an 
appropriate  punching  tool,  2.  Wash  off  remaining  electrolyte  in 
EMC,  3.  Drain  off  EMC  and  dry  the  anode  for  10  min  in  vacuum. 


gained  by  merging  electrochemical  impedance  spectroscopy  for  /  >  1  Hz  and  time 
domain  measurements  for/<  1  Hz. 
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As  the  LiFePCU/aluminum  sheets  are  always  coated  on  both 
sides  with  cathode  material,  one  coating  has  to  be  removed  by 
a  dissolving  liquid  and  a  cleaning  medium.  Hence,  the  preparation 
of  cathodes  comprises  this  additional  step.  This  step  also  necessi¬ 
tates  an  extensive  drying  for  at  least  3  h  at  100  °C  in  vacuum  in 
order  to  remove  the  entire  dissolving  liquid.  Having  been  prepared 
this  way,  anodes  and  cathodes  are  stored  under  glovebox  atmo¬ 
sphere  and  assembled  in  experimental  cells. 

2.2.  Experimental  cells  with  reference  electrode 


modeled  and  experimentally  verified  in  [25-27],  However,  the 
reference  electrode  potential  changes  slightly  at  small  currents. 

2.3.  Conducted  measurements 

All  measurements  were  carried  out  using  a  Solartron  1400E  cell 
test  system  with  Scribner  Multistat  software.  For  an  improved 
voltage  resolution  during  time  domain  measurements,  additionally 
an  Agilent  34970A  and  34901 A  were  used.  A  standard  temperature 
of  25  °C  was  adjusted  for  all  measurements  and  controlled  by 
a  Weiss  WK1  180  climate  test  chamber. 


The  experimental  cells  made  of  harvested  18650  cell  materials 
are  combined  in  commercially  available  EL-cell  housings  [20]  with 
a  connection  for  reference  electrodes  (ECC-Ref).  In  this  study,  all 
experimental  full  cells  (anode/cathode)  are  assembled  using 
a  reference  electrode  for  electrochemical  impedance  spectroscopy 
measurements  (Fig.  3),  whereas  all  half  cells  (anode/lithium  metal 
or  cathode/lithium  metal)  are  assembled  without  reference  elec¬ 
trode  for  time  domain  measurements.  The  half  cells  contain  one 
1.2  mm  thick  glass  fiber  separator  whereas  the  full  cells  contain 
four  220  pm  glass  fiber  separators,  two  on  each  side  of  the  reference 
electrode.  The  electrolyte  is  based  on  a  EC:EMC  1:1  solvent  mixture 
and  contains  a  1  M  mixture  of  LiPF6:LiC104  1 :9  as  conducting  salt. 
This  mixture  is  reported  to  prevent  Fe-dissolution  from  LiFePCU  and 
to  passivate  the  aluminum  current  collector  [21,22],  Exclusively  for 
time  domain  measurements,  5%  of  vinylene  carbonate  (VC)  is  added 
which  further  improves  the  cell  stability  for  the  extended 
measurement  time  [23,24], 

The  reference  electrode  used  in  this  study  is  a  V2A  steel  mesh 
(195  pm  mesh  size,  70  pm  thickness).  This  choice  was  made 
because  a  meshed  geometry  ensures  a  homogenous  current 
distribution  in  the  liquid  electrolyte.  This  arrangement  avoids 
inductive  artifacts  or  cross  talk  between  anode  and  cathode,  as 


LiFeP04 

:  electrolyte 
*■  separator 
reference 
electrode 
graphite 


Fig.  3.  Set-up  of  experimental  cell  for  impedance  measurements  via  reference 
electrode. 


2.3.1.  18650  Cell  test  procedure 

1.  Stability  of  the  18650  cell  was  ensured  by  10  times  cycling, 
following  a  CCCV  procedure  with  a  cutoff  voltage  of  3.6  V  and 
2.0  V,  respectively.  Subsequently,  the  state  of  charge  (SOC)  was 
varied  from  90%  to  0  %  SOC  in  10%  steps.  The  SOC  is  defined  by 
the  discharge  capacity  under  standard  discharge  conditions. 

2.  Time  domain  measurements:  the  18650  cell  was  equilibrated 
before  each  pulse  during  a  relaxation  phase  of  40  h.  Afterward, 
the  cell  was  excited  by  discharge  pulses  of  1C  height  and  10  s 
length.  The  relaxation  phase  for  the  following  40  h  was 
analyzed  by  Fourier  transformation  according  to  [2],  The 
measurement  procedure  for  one  SOC  is  shown  in  Fig.  4.  The 
subsequent  evaluation  of  impedance  (E1S)  and  time  domain 
measurements  is  described  in  the  next  chapter. 

3.  E1S  measurements:  potentiostatic  impedance  measurements 
were  conducted  for  the  high  and  medium  frequency  domain 
from  100  kHz  to  5  mHz.  The  voltage  amplitude  was  hereby 
chosen  to  be  10  mV  (100  kHz-10  Hz)  and  5  mV  (10  Hz-5  mHz). 
The  impedance  measurements  were  conducted  for  each  SOC 
after  the  time  domain  measurement. 

2.3.2.  Experimental  cell  test  procedure 

1.  Stability  of  experimental  cells:  the  assembled  experimental 
cells  were  cycled  for  20  cycles  following  a  CCCV  protocol  in 
order  to  conduct  the  formation  process.  The  cutoff  voltage  for 
full  cells  was  the  same  as  for  the  18650  cells,  namely  3.6  V  and 
2.0  V.  For  anode  half  cells  it  was  0.01  V  and  0.8  V  whereas  the 
cathode  half  cells  cycled  between  2.7  V  and  3.7  V.  Afterward, 
impedance  and  time  domain  measurements  were  recorded  for 
varying  SOC.  The  impedance  measurements  were  conducted 
for  full  cells  with  reference  electrode  whereas  the  time  domain 
measurements  were  conducted  for  full  cells  without  reference 
electrode  and  anode  and  cathode  half  cells  with  lithium 
counter  electrode. 
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2.  Time  domain  measurements  (half  cells  and  full  cells  without 
reference):  time  domain  measurements  were  conducted  for 
80%  SOC.  Fig.  4  shows  the  measurement  sequence  for  a  time 
domain  measurement.  For  experimental  cells,  the  discharge 
pulse  was  also  chosen  as  1C/10  s  pulse  and  the  voltage  relax¬ 
ation  phase  was  40  h  as  well. 

3.  EIS  measurements  (full  cells  with  reference):  potentiostatic 
impedance  measurements  were  conducted  for  the  high  and 
medium  frequency  range  from  100  kHz  to  10  mHz.  The  voltage 
amplitude  was  hereby  chosen  to  be  10  mV  (100  kHz-100  Hz) 
and  5  mV  (100  Hz-10  mHz).  The  SOC  was  varied  from  100% 
to  0%  SOC  in  10%  SOC  steps. 

3.  Analysis  of  EIS  and  TDM  measurement  results 

3.1.  Electrochemical  impedance  spectroscopy  (EIS)  measurements 


i  i 

“ 

■  experimental  full  < 

■  cathode  via  refer* 

■  anode  via  referen 

cell 

ce 

Fig.  6.  Real  part  and  imaginary  part  of  Kramers  Kronig  resii 
cell  (black)  and  anode  (red)  and  cathode  (blue)  reference  electrode  measurements.  (For 
interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader  is  referred  to 
the  web  version  of  this  article.) 


3.1.1.  Measurement  data  quality  (Kramers  Kronig  residuals) 

At  first,  the  data  quality  of  our  impedance  measurements  is 
assessed  by  calculating  the  real  and  the  imaginary  part  of  the 
Kramers  Kronig  (KK)  residuals  [28],  Kramers  Kronig  residuals  of 
a  full  cell  impedance,  measured  in  an  experimental  cell,  are 
compared  to  those  of  the  18650  cell.  As  Fig.  5  demonstrates,  the 
residuals  of  both  measurements  are  below  0.3%  or  even  0.2%  within 
the  frequency  range  from  0.1  Hz  to  5  kHz.  This  data  quality  is 
mandatory  for  a  subsequent  DRT  evaluation.  However,  the  Kramers 
Kronig  residuals  increase  to  1%  for  high  frequencies  /  >  10  kHz,  as 
the  measurement  range  is  changed  automatically  in  the  Solartron 
measurement  equipment.  This  causes  a  small  misalignment  in  the 
impedance  spectrum  which  cannot  be  described  by  the  Kramers 
Kronig  transformation.  Furthermore,  the  inductive  tail  shows  no 
ideal  inductive  behavior  which  also  causes  increasing  residuals. 

The  Kramers  Kronig  test  is  of  the  greatest  importance  for 
impedance  measurements  between  a  working  electrode  and 
a  reference  electrode,  as  this  potential  is  not  necessarily  stable  over 
time.  This  instability  transgresses  the  EIS  pre-condition  of  time 
invariance  and  limits  the  validity  of  measurement  results.  The 
influence  of  potential  change  is  here  evaluated  by  the  size  of 
Kramers  Kronig  residuals  which  is  demonstrated  in  Fig.  6.  Herein, 
the  KK  residuals  for  the  full  cell  impedance  are  plotted  together 
with  the  KK  residuals  for  the  half  cell  impedances,  i.e.,  anode  and 
cathode  versus  reference  electrode  measurements.  As  before,  KK 
residuals  at  high  frequencies  (f  >  10  kHz)  deviate  to  higher  values 
caused  by  the  automatically  changing  measurement  range.  The  KK 
residuals  at  medium  frequencies  (1  Hz  >/>  5  kHz)  are  well  below 
0.3%,  thus  proving  a  very  good  data  quality.  However,  the  KK 
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Fig.  5.  Real  part  and  imaginary  part  of  Kramers  Kronig  residuals  for  experimental  full 
cell  (black)  and  18650  cell  (blue)  measurements.  (For  interpretation  of  the  references 
to  colour  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


residuals  for  measurements  via  a  reference  electrode  increase 
significantly  to  1%  or  more  for  low  frequencies  (f  <  1  Hz),  which  is 
attributed  to  the  above  mentioned  instabilities  over  time.  Please 
recall  from  chapter  2.3,  that  pulse  measurements  were  applied  on 
half  cells  without  reference  electrode  in  order  to  avoid  the  influ¬ 
ence  of  shifting  reference  electrode  potential. 

In  this  study,  therefore,  EIS  measurements  are  evaluated  only  in 
the  frequency  range  100  kHz  >/>  1  Hz.  Below  these  frequencies, 
impedance  data  are  evaluated  only  from  time  domain 
measurements. 

3.1.2.  EIS  data  evaluation  by  the  distribution  of  relation  times  (DRT) 

In  general,  EIS  is  especially  useful  if  the  system  performance  is 

governed  by  coupled  processes  each  proceeding  at  a  different  rate. 
However,  their  number  and  physical  nature  are  not  always  identi¬ 
fied  unambiguously.  For  electrodes  with  a  complex  microstructure, 
several  reaction  steps  contribute  which  are  hardly  separable  in  the 
impedance  response,  as  their  relaxation  frequencies  widen  out  and 
overlap  in  the  impedance  spectrum. 

This  problem  can  be  avoided,  e.g.,  by  the  following  alternative 
approach  for  analyzing  impedance  spectra.  This  procedure  was 
developed  for  the  impedance  assessment  of  solid  oxide  fuel  cells 
and  was  presented  in  detail  in  [4-6].  Recently,  it  has  been  applied 
for  the  first  time  for  the  analysis  of  experimental  lithium  ion  cells 
[7—9],  The  basic  concept  of  DRT  evaluation  is  described  in  the  next 
subchapter. 

3.1.3.  Theory  of  distribution  function  of  relaxation  times  (DRT) 

The  DRT-method  uses  the  fact,  that  every  impedance  function 
that  obeys  the  Kramers— Kronig  relations  can  be  represented  as  an 
infinite  number  of  infinitesimal  small  differential  RC-elements  [3], 
This  is  common  in  system  theory  and  also  holds  true  for  electro¬ 
chemical  systems  like  lithium  ion  batteries. 

For  simplification,  Fig.  7a  shows  the  DRT  for  two  ideal  processes 
which  can  be  represented  by  two  RC  elements.  The  RC  elements 
obey  the  following  equation: 

7  7  7  R|  f?2 

Complete  =  4?C1  +4tC2  =  )  +jll)R]C]  +  )  +j(l)R2C2 

e  _ ^ _ i--  -  % _  (1) 

f  +JWT!  1  +  jm2 

There  are  two  parameters  which  characterize  each  RC  element 
and  therefore  the  underlying  physical  process.  The  first  one  is  R, 
giving  the  polarization  of  the  corresponding  process.  The  second 
one  is  the  time  constant  t,  representing  the  relaxation  frequency 
or  rather  relaxation  time  of  the  process.  It  is  equal  to 
fr  =  RiCi  =  (1/271-/0. 
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Fig.  7.  DRT  for  a)  ideal  processes,  represented  by  RC-elements  and  b)  real  processes, 
represented  by  RQ-elements. 


obtained  by  division  of  the  Fourier  coefficients  [2],  In  this  study,  we 
use  a  current  pulse  as  excitation  signal  and  analyze  the  system 
response  during  and  after  the  excitation. 

An  ideal  current  pulse  stimulates  low  frequencies  and  is  there¬ 
fore  adequate  to  determine  an  impedance  response  for  several 
sufficiently  low  frequencies  by  one  single  measurement.  The 
frequency  range  of  a  time  domain  measurement  is  only  restricted 
by  the  limited  measurement  time  for  low  frequencies  and  the 
limited  sample  rate  for  high  frequencies. 

The  differential  capacity  AC  complicates  the  DRT  evaluation  of 
a  data  set  assembled  via  EIS  [8],  However,  in  TDM  it  is  easily 
determined  by  the  charge  of  the  current  pulse  and  the  voltage  gap 
A U  before  and  after  the  current  pulse  AQ,  It  is  defined  as 
AC  =  (AQ/AU)  and  can  be  directly  subtracted  from  the  impedance 
spectrum  measured  via  TDM. 

4.  Results  and  discussion 


In  the  given  case,  the  process  parameters  R  and  t  can  be  directly 
read  from  the  DRT  by  height  and  location  of  the  two  Dirac  peaks.  In 
the  following  figures,  we  use  the  relaxation  frequency/  =  (1  / 2m ) 
instead  of  the  relaxation  time  t  for  the  x-values. 

Advanced  electrodes  with  a  complex  microstructure  do  not  have 
one  time  constant  but  a  distribution  of  time  constants.  An  example 
for  such  a  real  system  is  given  in  Fig.  7b.  Each  polarization  process  is 
now  characterized  in  the  DRT  by  the  center  frequency  of  its  peak 
which  is  the  corresponding  relaxation  frequency  and  by  the  shape 
of  its  peak.  The  area  under  the  peak  enclosed  to  the  peak  frequency 
corresponds  to  the  polarization  of  the  loss  process.  The  peak  shape 
of  each  peak  depends  on  the  nature  of  the  underlying  process. 
Some  DRT  examples  for  commonly  used  impedance  elements  as 
Warburg  element  are  shown  in  [5], 

3.1.4.  DRT  calculation 

The  relation  between  DRT  and  impedance  curve  is  given  [4]  by 

m  =  Ro  +  Zpo,(<o)  =R0  +  RPo,  ■  J  (2) 

o 

The  term  (g(r)/l  +jcn)dx  specifies  the  fraction  of  the  overall 
polarization  with  relaxation  times  between  t  and  t  +  dr.  Rq  is  the 
ohmic  resistance  and  Rp0i  the  overall  polarization.  From  linear 
system  theory  it  is  well  known  that  the  impedance  of  each  entirely 
capacitive  electrical  system  can  be  transformed  into  the  form  in  Eq. 
(2).  However,  for  a  numerical  DRT  calculation,  only  a  finite  number 
of  RC-elements  can  be  used.  This  finite  number  of  RC-elements  can 
be  represented  by 

W">  -  (3) 

giving  the  overall  polarization  resistance. 

The  mathematical  problem  with  the  DRT  approach  arises  from 
the  inversion  of  Eq.  (3)  which  is  necessary  in  order  to  extract  g  (t) 
from  the  measured  impedance  data  Zpo\  (w).  This  problem  is  known 
to  be  ill-posed  and  requires  special  methods  to  be  solved  in  order  to 
avoid  false  peaks  and  oscillations  [4,29,30], 

3.2.  Time  domain  measurements  (TDM) 

For  time  domain  measurements,  a  pulse  signal  and  the  corre¬ 
sponding  system  response  is  transformed  via  Fourier  trans¬ 
formation  to  the  frequency  domain,  and  the  impedance  spectrum  is 


4.1.  Impedance  spectra  measured  by  EIS 

4.1.1.  18650  Cell 

Fig.  8a  shows  the  impedance  spectra,  measured  by  EIS  only,  of 
the  examined  18650  cell  for  90%  SOC.  The  impedance  is  scaled  to 
the  electrode  area  (Q  cm2),  which  enables  a  direct  comparison  to 
the  impedance  spectra  of  experimental  cells  presented  in  chapter 

4.1.2. 

A  capacitive  diffusion  branch  for  low  frequencies,  one  flat 
semicircle  for  medium  frequencies  and  the  typical  inductive 
behavior  for  high  frequencies  are  immediately  observable  from 
the  impedance  spectrum  in  Fig.  8a.  However,  the  transformation 
of  these  data  into  the  space  of  relaxation  times,  shown  in  Fig.  8b, 
opens  up  to  details  especially  in  the  medium  frequency  range. 
Several  peaks  appear  between  1  Hz  and  5  kHz  which  all  add  up 
to  one  single  semicircle  in  the  corresponding  Nyquist  plot. 
However,  frequencies  <1  Hz  cannot  be  de-convoluted  by  DRT 
without  using  a  pre-processing  step,  as  the  DRT  calculation 
requires  a  convergent  impedance  curve.  Therefore,  the  subtrac¬ 
tion  of  the  capacitive  diffusion  branch  becomes  indispensable. 
The  feasibility  of  this  approach  was  worked  out  in  detail  in  our 
recent  paper  [8], 
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Fig.  9.  a)  Impedance  spectra  and  b)  DRT  of  18650  cell  for  SOC-variation  (90%-0%). 


In  this  study,  we  present  a  new  approach,  as  the  low  frequency 
range  will  now  be  investigated  by  time  domain  measurements, 
which  works  without  any  pre-processing.  For  high  frequencies,  the 
evaluation  is  limited  by  the  inductive  characteristic  of  the 
measured  system.  The  latter  overlaps  the  high  frequency  loss 
processes  which  might  occur  in  the  impedance  response  and 
prevent  an  interpretation  using  the  DRT  for  frequencies  above 
5  kHz.  These  frequencies  are  therefore  only  analyzed  in  the 
experimental  cells,  as  their  inductance  is  much  lower. 

Furthermore,  the  18650  cell  is  characterized  by  varying  SOC 
from  90%  to  0%,  as  shown  in  Fig.  9a.  Obviously,  its  impedance 
changes  not  much  with  SOC,  as  only  for  low  SOC  values  the 
polarization  increases  slightly.  The  same  conclusion  holds  for  the 
calculated  DRTs,  as  depicted  in  Fig.  9b. 


4.1.2.  Experimental  cells  with  reference  electrode 

Fig.  10  compares  the  impedance  spectra  and  DRTs  for  the 
experimental  full  cell  at  100%  SOC  with  the  impedance  measure¬ 
ment  and  DRT  evaluation  via  reference  electrode.  The  full  cell 
impedance  spectrum  in  Fig.  10a  reveals,  beside  the  capacitive 
diffusion  branch  for  low  frequencies  and  the  inductive  behavior  for 
frequencies  above  100  kHz,  two  characteristic  semicircles. 
Regarding  the  impedance  curves  of  anode  and  cathode  via  refer¬ 
ence  electrode  in  Fig.  10b /c,  two  semicircles  can  be  assigned  to  the 
anodic  and  one  semicircle  to  the  cathodic  electrode.  The  ohmic 
resistance  splits  up  equally  for  both  electrodes,  which  indicates  an 
optimal  position  of  the  reference  electrode.  Further,  no  inductive 
loops  or  other  typical  artifacts  can  be  observed  for  the  impedance 
measurements  via  reference  electrode  [25,26]. 

Fig.  10  d,  e  and  f  show  the  calculated  DRTs  in  order  to  give 
a  better  separation  of  anode  and  cathode  contributions.  The 
measured  frequencies  are  divided  into  three  characteristic 
frequency  ranges,  denominated  as  I,  II,  and  III.  Frequency  range  1 
includes  the  high  frequencies  between  5  kHz  and  100  kHz.  Herein, 
two  characteristic  peaks  are  clearly  distinguishable  in  the  calcu¬ 
lated  DRT  of  the  full  cell.  The  DRTs  of  anode  and  cathode  are  rather 
similar  within  this  frequency  range;  therefore  both  electrodes 
contribute  to  the  impedance  equally.  Frequency  range  II  covers 
frequencies  between  10  Hz  and  5  kHz.  In  the  DRT,  the  curve  shape 
of  the  full  cell  indicates  the  existence  of  up  to  four  characteristic 
time  constants,  which  is  a  first  hint  of  either  four  individual  loss 
processes  or  less,  if  a  loss  process  with  several  coupled  time 
constants  takes  place  here.  In  this  case,  an  inspection  of  the  refer¬ 
ence  electrode  measurements  is  supportive  again,  clearly  indi¬ 
cating  that  the  anode  dominates  the  impedance  contributions  in 
the  entire  frequency  range  II. 

Frequency  range  III,  assigned  to  frequencies  below  10  Hz,  is  only 
to  be  analyzed  by  time  domain  measurements  in  Section  4.2.  This 
procedure  is  required,  because  EIS  impedance  data  contains 
a  differential  capacity,  which  will  cause  misleading  peaks  in  the 
DRT  calculation.  We  have  already  evaluated  this  in  detail  and 


Fig.  10.  Left:  Impedance  spectra  measured  by  EIS  of  a)  an  experimental  full  cell  (black),  b)  cathode  (blue)  and  c)  anode  (red)  determined  via  reference  electrode;  right:  d),  e)  an 
calculated  corresponding  DRTs,  all  at  100%  SOC.  (For  interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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introduced  a  pre-processing  procedure,  as  shown  in  [8],  to  over¬ 
come  this  problem.  Alternatively,  as  stated  in  Section  3.2,  this 
frequency  range  III  can  be  analyzed  directly  by  time  domain 
measurements. 

Now  we  have  a  look  at  the  SOC  variations  between  100%  and  0% 
at  25  °C,  which  are  displayed  only  as  DRTs  in  Fig.  11.  First,  the  DRTs 
in  range  II  are  quite  comparable  between  18650  cell  (Fig.  9)  and 
experimental  cell.  There  are  some  smaller  variations  in  the  shape  of 
the  calculated  DRTs  at  SOC  >10%  within  range  II,  whereas  range  I 
remains  constant.  Noteworthy  in  range  II  is  the  shape  of  the  DRT  at 
SOC  =  0%,  which  corresponds  to  an  open  circuit  potential  of  2.6  V. 
Here,  the  predominant  peak  of  the  full  cell  increases  significantly  in 
height,  changes  its  shape  and  its  maximum  shifts  to  a  more  than 
one  order  of  magnitude  lower  frequency.  The  measurements  via 
reference  electrode  in  Fig.  llb/c  reveal  that  this  effect  is  caused 
exclusively  by  the  anode,  whereas  the  cathode  remains  unaffected. 
Please  note,  that  we  shifted  the  lower  border  of  frequency  range  II 
from  10  Hz  in  Fig.  10  to  1  Hz  in  Fig.  11  in  order  to  include  the  SOC 
dependent  peak. 

4.1.3.  Comparison  of  experimental  cell  and  18650  cell  results  (EIS 
and  DRT) 

Fig.  12a  compares  EIS  measurements  of  (i)  the  18650  cell  (blue 
solid  line)  at  SOC  =  0%  (3.08  V)  with  (ii)  the  experimental  cell  (black 
solid  line)  at  SOC  =  10%  (2.98  V)  and,  (iii)  the  same  spectrum  of  the 
18650  cell  (blue  dotted  line),  but  shifted  toward  lower  frequencies 
that  the  capacitive  branches  of  both  cells  overlap.  The  corre¬ 
sponding  DRTs  of  both  cell  types  are  given  in  Fig.  12b.  A  different 
SOC  for  both  cell  types  was  chosen,  because  the  DRT  of  the 
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Fig.  11.  DRTs  calculated  from  EIS  measurements  of  a)  experimental  full  cell,  b)  cathode 
and  c)  anode  measured  via  reference  electrode  at  SOC-variation. 


experimental  cell  steeply  increases  at  SOC  =  0%  within  frequency 
range  II  (see  Fig.  11),  and  because  open  circuit  potentials  among 
both  cells  should  be  as  close  as  possible. 

First  of  all,  the  impedance  polarization  of  the  experimental  cell 
measurement  is  in  the  order  of  the  18650  cell,  when  both  curves  are 
normalized  on  the  electrode  area.  However,  the  shape  of  the 
impedance  spectra  differs  noticeably,  as  for  the  experimental  cell  (i) 
the  resistance  R0  is  higher,  and  (ii)  the  polarization  resistance  Rpoi 
shows  two  arcs  for  medium  frequencies.  These  characteristics 
evolve  from  the  difference  in  the  inductivity  of  both  type  of  cells. 
This  is  concluded  from  the  x-axis  interception  frequency,  which  is 
at  3  kHz  for  the  18650  cell  and  at  100  kHz  for  the  experimental  cell. 
Therefore,  Fig.  12a  also  shows  a  shift  of  the  18650  cell  impedance 
curve  toward  higher  real  part  values,  until  the  capacitive  branch  of 
both  cell  types  are  coincident.  This  demonstrates  that  the  high 
frequency  part  of  the  18650  cell  impedance  is  governed  by  the 
inductivity  of  current  collectors,  which  superimposes  the  polari¬ 
zation  losses  taking  place  at  the  same  frequencies.  More  informa¬ 
tion  is  gained  from  the  comparison  of  both  calculated  DRTs, 
displayed  in  Fig.  12b. 

Range  I  is  only  accessible  by  experimental  cells  with  a  low 
inductivity.  As  already  explained  in  Section  4.1.2,  two  characteristic 
peaks  are  clearly  distinguishable,  and  anode  and  cathode 
contribute  to  the  impedance  equally.  It  is  deduced,  that  these  high 
frequency  contributions  occur  as  well  in  the  18650  cell,  but  they  are 
included  in  the  ohmic  resistance  Ro  of  the  18650  cell.  The  DRTs  in 
range  II,  with  frequencies  correlating  to  the  second  semicircle  of  the 
impedance  spectrum  of  the  experimental  cell,  show  very  good 
agreement  in  the  DRTs.  Both  type  of  cells  show  four  peaks  located 
at  the  same  characteristic  frequencies.  Merely  one  peak  shows 
different  size  of  polarization.  This  might  result  from  the  electrode 
preparation  procedure  or  the  use  of  a  different  electrolyte.  As 
already  discussed  in  Section  4.1.2,  the  anode  dominates  the 
impedance  response  over  the  entire  frequency  range. 

Frequency  range  III,  assigned  to  frequencies  lower  than  10  Hz,  is 
interpreted  only  by  time  domain  measurements  in  Section  4.2. 


Fig.  12.  Comparison  of  a)  EIS  measurements  of  experimental  full  cell  (black)  and  18650 
cell  (blue)  and  b)  corresponding  calculated  DRTs,  both  displayed  for  the  18650  cell  at 
SOC  =  0 %  (3.08  V)  and  for  the  experimental  cell  at  SOC  =  10%  (2.98  V).  (For  inter¬ 
pretation  of  the  references  to  colour  in  this  figure  legend,  the  reader  is  referred  to  the 
web  version  of  this  article.) 
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4.2.  Impedance  spectra  measured  by  TDM 

Impedance  spectra  derived  from  time  domain  measurements 
and  the  correlated  DRTs  shown  in  the  following  chapter  are  scaled 
to  the  electrode  area,  which  facilitates  a  comparison  of  the  18650 
cell  with  the  experimental  cell.  TDM  measurements  cover  the 
frequencies  between  2  pHz  and  100  mHz. 

4.2.1.  18650  Cell 

Fig.  13  shows  the  current  pulse  (lower  part)  with  a  duration  of 
10  s,  which  is  put  on  both  cell  types  as  excitation  signal  (with 
different  current  rate)  and  the  voltage  response  for  the  18650  cell 
(upper  part)  for  a  time  domain  measurement  at  80%  SOC.  The  cell 
voltage  decreases  during  the  discharge  current  pulse.  A  small  part 
of  the  relaxation  phase  is  enlarged  in  order  to  show  the  most 
significant  voltage  response.  Fig.  14  delivers  the  impedance  spec¬ 
trum  as  well  as  the  DRT  calculated  from  the  time  domain  response. 

The  impedance  spectrum  in  this  frequency  range  is  composed  of 
one  large  semicircle,  which  is  much  larger  than  those  presented  in 
the  previous  chapter.  The  DRT  shows  a  decreasing  peak  sequence, 
which  is  typical  for  a  polarization  loss  describable  by  a  Warburg 
element  [5],  Furthermore,  the  DRT  reveals  one  additional  large 
peak  at  very  low  frequencies  off  =2  pHz  which  superimposes  the 
Warburg  impedance. 


a) 


Fig.  14.  a)  Impedance  spectrum  measured  by  TDM  and  b)  corresponding  calculated 
DRT  of  an  18650  cell  at  80%  SOC. 


4.2.2.  Experimental  cells  without  reference  electrode 

Fig.  15  shows  the  impedance  spectra  and  the  corresponding 
DRTs  of  the  experimental  full  cell  and  the  two  half  cell  measure¬ 
ments  at  80%  SOC  and  25  °C.  The  impedance  spectra  of  full  cell  and 
cathode,  as  for  the  18650  cell,  are  composed  of  one  large  semicircle, 
whereas  the  anode  displays  a  negligible  small  contribution  in  this 
frequency  range.  The  DRTs  of  full  cell  and  cathode  show 
a  decreasing  peak  sequence,  describable  by  a  Warburg  element.  It  is 
concluded,  that  the  cathode  dominates  the  impedance  contribu¬ 
tions  in  this  frequency  range. 

4.2.3.  Comparison  of  experimental  cell  and  18650  cell  results  (TDM) 

Now  we  compare  the  results  compiled  by  TDM  and  DRT  analysis 

in  Fig.  16.  The  impedance  spectra  of  the  18650  cell  and  the  experi¬ 
mental  cell  differ  in  polarization  size  and  characteristic  shape,  which 
holds  for  the  corresponding  DRTs  as  well.  Herein,  the  measured 
frequency  range  is  now  divided  in  two,  range  III  and  range  IV.  Range 
III  was  already  previously  assigned  to  the  frequencies  between  5  pHz 
and  10  Hz.  Herein,  the  DRTs  of  experimental  cell  and  18650  cell  are 
appropriate  in  polarization  magnitude  and  characteristic 
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Fig.  13.  First  part  of  current  step  and  voltage  response  for  the  time  domain 
measurement  of  an  18650  cell  for  80%  SOC. 
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frequencies.  As  already  stated  in  the  previous  section,  these 
contributions  are  solely  caused  by  the  cathode. 

For  the  newly  added  range  IV,  only  the  DRT  of  the  18650  cell 
shows  an  additional  peak  near  2  pHz.  Thereby,  the  DRT  delivers 
very  helpful  information,  as  this  peak  explains  the  difference 
among  both  impedance  spectra. 

4.3.  Physical  interpretation 

As  already  indicated  in  Sections  4.1.3  and  4.2.3,  the  impedances 
measured  by  E1S  and  TDM  for  an  experimental  cell  are  supportive 
for  the  physical  interpretation  of  the  impedance  of  an  18650  cell. 

4.3.1.  Impedance  spectra  measured  by  EIS 

The  resistance  value  Rq  of  the  18650  cell,  determined  in  the  high 
frequency  range  above  5  kHz  (range  I),  is  composed  of  electrode 
polarization  processes  and  the  ohmic  resistance  of  the  electrolyte 
(Fig.  12).  However,  the  origin  of  the  high  frequency  processes  can  be 
identified  easily  in  experimental  cells  (Fig.  10).  We  assume,  that  the 
polarization  losses  at  high  frequencies  originate  from  the  two 
individual  contact  resistances  between  anode/Cu  and  cathode/Al. 
This  physical  interpretation  was  determined  in  previous  investi¬ 
gations  for  a  LiFePC>4-cathode/Al  contact  resistance  [8,31]  and  is 
here  assumed  to  be  transferable  to  a  graphite-anode/Cu  contact 
resistance.  Therefore,  the  resistance  value  R0  of  the  18650  is 
attributed  to  (i)  the  two  contact  resistances  between  both  elec¬ 
trodes  and  the  respective  current  collector  and  (ii)  the  ohmic 
resistance  of  the  electrolyte. 

In  the  medium  frequency  range  (range  II),  the  impedance  curves 
for  experimental  cell  and  18650  cell  are  very  similar  (Fig.  12).  This 
fact  allows  transferring  the  conclusions  of  the  measurements  via 
reference  electrodes  to  the  18650  cell.  Therefore,  the  anode  domi¬ 
nates  the  impedance  response  of  the  18650  cell  over  several 
frequency  decades  (10  Hz— 5  kHz,  Fig.  10).  As  the  DRT  reveals  anode 
peaks  stretched  over  2.5  decades,  this  indicates  a  mixed  contribu¬ 
tion  consisting  of  SEI  in  combination  with  charge  transfer.  Herein, 
the  SOC-variation  indicates  a  strong  dependency  especially  at  the 
lower  SOC  values  for  the  DRT  peak  with  the  slowest  time  constant 
(Fig.  11 ).  This  DRT  peak  is  therefore  assigned  to  the  charge  transfer 
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Fig.  15.  Left:  Impedance  spectra  measured  by  TDM  and  right:  corresponding  DRTs  of  (a,  b)  experimental  full  cell  (black),  (c,  d)  cathode  (blue)  and  (e,  f)  anode  (red)  half  cells  at  80% 
SOC.  (For  interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


process  anode/electrolyte,  whereas  the  faster,  SOC  independent 
peaks  are  assigned  to  the  SEI  layer. 

4.3.2.  Impedance  spectra  measured  by  TDM 

The  time  domain  measurements  of  experimental  cells  (Fig.  15) 
point  to  a  Warburg  type  impedance,  which  delivers  a  decreasing 
peak  sequence  in  the  DRT  representation  [5],  This  characteristic 


f/Hz 

Fig.  16.  Comparison  of  experimental  cell  (black)  and  18650  cell  (blue)  via  a)  imped¬ 
ance  spectra  measured  by  TDM  and  b)  corresponding  calculated  DRTs  at  80%  SOC.  (For 
interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader  is  referred  to 
the  web  version  of  this  article.) 


points  to  a  solid  state  diffusion  taking  place  inside  the  active 
material.  As  shown  for  the  experimental  cell,  the  cathode  is  most 
relevant  for  this  polarization  contribution.  Therefore,  the  solid  state 
diffusion  taking  place  in  the  LiFePCU  cathode  of  the  18650  cell  is  the 
dominant  polarization  loss  over  the  entire  frequency  range  III. 

Below  5  pHz  (range  IV),  the  DRT  reveals  an  additional  peak 
occurring  only  in  the  18650  cell,  but  not  present  in  the  experi¬ 
mental  cell  (Fig.  16).  Hence,  this  polarization  process  remains 
unassignable  to  anode  or  cathode,  as  this  DRT  peak  (a)  does  not 
exist  in  the  experimental  cell  or  (b)  cannot  be  measured  due  to 
stability  problems  at  these  very  low  frequencies. 

Assuming  cause  (a),  the  physical  origin  may  be  related  with 
a  homogenization  process  inside  the  entire  volume  of  the  electrode 
layers.  Probably,  we  have  uncovered  lithium  concentration  gradi¬ 
ents  among  the  electrodes  of  the  18650  cell.  Naturally,  such 
concentration  gradients  do  not  exist  in  the  experimental  cells 
having  only  a  very  limited  electrode  area. 

5.  Conclusions 

The  physical  origin  of  the  impedance  of  a  commercial  18650  cell 
was  investigated  within  a  wide  frequency  range  (100  kHz— 2  pHz) 
by  combining  electrochemical  impedance  spectroscopy  (EIS)  and 
time  domain  measurements  (TDM).  In  parallel,  the  original  18650 
cell  was  opened  and  impedance  response  of  anode  and  cathode  was 
separated  in  an  experimental  test  set-up  using  a  reference  elec¬ 
trode.  Moreover,  the  different  kinds  of  impedance  data  sets  were 
de-convoluted  in  the  space  of  relaxation  times  (DRT).  This  multistep 
approach,  which  is  not  tied  to  a  particular  electrochemistry, 
allowed  the  identification  of  the  dominating  physical  processes 
being  hidden  in  the  impedance  spectrum  of  the  original  18650  cell: 

-  range  I  from  100  kHz  to  5  kHz  (measured  by  EIS):  inductive 
behavior,  most  probably  originating  from  the  internal  current 
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collectors,  dominates  the  impedance  response  of  the  18650  cell 
in  this  frequency  range.  Experimental  cell  measurements  dis¬ 
closed  that  the  electrolyte  contribution  to  the  ohmic  resistance 
of  the  18650  cell  is  overestimated,  as  not  only  the  electrolyte  as 
usually  assumed  [1],  but  also  polarization  losses  originating 
from  contact  resistances  anode/Cu  and  cathode/Al  contribute 
to  the  “ohmic”  resistance,  and  do  not  change  with  SOC. 

-  range  II  from  5  kHz  to  10  Hz  (measured  by  E1S):  the  anode 
dominates  the  impedance  response  of  the  18650  cell  over  the 
entire  frequency  range.  As  the  DRT  de-convolutes  several 
anode  peaks  with  time  constants  stretched  over  2.5  decades, 
this  indicates  a  mixed  contribution  of  the  interface  SEI/anode 
in  combination  with  a  charge  transfer  reaction  to  the  electro¬ 
lyte.  The  low  frequency  peak  assigned  to  charge  transfer 
increases  strongly  with  decreasing  SOC. 

-  range  III  from  10  Hz  to  5  pHz  (measured  by  TDM):  the  solid 
state  diffusion  taking  place  in  the  LiFeP04  cathode  of  the  18650 
cell  is  the  dominant  polarization  loss  over  the  entire  frequency 
range.  This  is  an  unexpected  result,  however,  the  anode 
contribution  to  the  polarization  loss  seems  to  be  negligible 
small.  The  portions  of  cathode  and  anode  show  no  major 
changes  with  SOC. 

-  range  IV  at  /  <  5  pHz  (measured  by  TDM):  probably,  the 
impedance  response  of  the  18650  cell  matching  with  such 
a  time  constant  may  be  related  with  a  homogenization  process. 
It  is  possible  that  lithium  concentration  gradients  among  the 
entire  volume  of  the  electrodes  of  the  18650  cell  are  taking  place 
here.  This  frequency  range  was  only  analyzed  for  80%  SOC. 
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